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AB3TEACT 


Woven-soreen  or  orossed-rod  types  of  matrices 
may  be  of  interest  as  fuel— element  geometries  for 
oertaia  types  of  nuolear  reactors  as  well  as  in  many 
other  heat- transfer  and  mass-transfer  systems*  Basic 
heat-transfer  and  flow-friction  design  data  for  those 
geometries  are  presented  graphloally,  in  nondlmensional 
form,  for  a  range  of  Beynolds  numbers  from  5  to  100,000 
and  matrix  porosities  from  0*60  to  0»83.  Algebraic 
equations  which  adequately  represent  these  data  are 
also  given  cmd  these  may  be  used  for  Interpolation  and, 
at  least  tentatively,  for  a  limited  extrapolation  be¬ 
yond  the  range  of  porosities  covered  in  this  program. 
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OF  WOm-SOREEN  AND  OROSSEIW-'ROD  MATRIOBS 
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The  following  nomenolature  is  used  in  the  paper 

-  heat- transfer  area^  sq,  ft 

-  matrix  average  flow  area^  »  ^^frontai® 

“  matrix  frontal  area,  sq  ft 

-  speoifio  heat  at  constant  pressure^  Btu/(lb  deg  F) 

-  thermal  oapaoity  of  copper  cylinder,  Btn/deg  F 

-  diameter  of  rod  and  wire,  in®,  ft 

-  base  of  natural  logarithms,  2o?l8ooo 

o  factor  defined  in  Equation  (14) 

reoiprooal  of^proportionallty  factor  in  Newton's  second  law,  32,2 

(lb/)!)  («/••«*) 

-  average  mass  velooity  In  matrix,  based  on  lb/ (hr  sq  ft) 

-  maximum  mass  velocity  In  matrix,  based  on  minimum  free-flow  area, 
^^frontal® 

•»  unit  oonduotanoe  for  thermal  oonveotlon  heat  transfer,  Btu/(hr  sq 
ft  deg  F) 

-  \2nlt  thermal  conductivity,  Btu/(hr  sq  ft  deg  F) 

-  length  of  matrix  In  the  flow  direction,  ft 

-  porosity  of  matrix,  (vol  of  voldo)/(vol  of  matrix) 

-  pressure,  psf,  or  Inches  of  E^O 

-  pressure  differential,  psf,  or  inches  of  H^O 

-  hydraulic  radius.  Equations  (l)  and  (4)b  ^t 

-  mass-flow  rate,  Ibs/hr 
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one  TAB  □ 

Uriannounced  □ 

Justification 


By 

■ 

Oist.  ibutfon  / 

Availability  Codes 

Oist 

Avail  and/or 

Special 

-heat- transfer  area  per  unit  volume  of  matrixj  Equation  (3)»  ft 
6  -time^  hr}  and  angle  of  inolinatiom 
|X  -visoosityy  lib/(hr  ft) 

P  -mass  density;  pof 

^  -ratio  of  minimum  free-flow  to  frontal  area  ("fractional  opening  in 
table  l)»  dimensionless 

-wall  friction  drag  force  per  unit  of  wall  surface;  psf 
^  -screen  thickness pg  -  2d  x  ,  in*;  ft 

Dimensionless  Groups 
f  -Panning  friction  factor 

Cjj  -drag  coefficient  per  screen 

Np^  -prandtl  number;  Np^ 

Hp  -Heynolds  number  for  tube-llke  flow  or  interior  flow; 

E'p  -modified  Heynolds  number;  N*p"  *^H 

Hp(d)  -Reynolds  nvimber;  for  flow  over  a  body;  or  exterior  flow;  ^xaax*'' 

-Stanton  number;  ■  ^/®®p 

—longitudinal  pitch;  ratio  of  longitudinal  spaoiag  to  wire  or  rod 
diameter 

X  -transverse  pitch;  the  ratio  of  transverse  spacing  to  wire  or  rod 
^  diameter 


The  symbol  #  denotes  pounds  force  in  distinotion  to  lb  for  pounds 


mass. 


•denotes  defining  equation 


UTEODUOTIOI 


Matrices  of  the  woven-screeny  orossed-rod,  or  sphere-bed  types  may 
prove  to  be  of  Interest  as  fuel  element  geometries  for  gas-oooled,  diphenyl- 
oooled,  and  nonboiling  water  reaotors«  Fine  structure  matrloes  offer  the 
advantages  of  (a)  high  transfer  area  densities,  sq  ft  of  surface/  ou  ft  of 
volume)  (b)  excellent  heat- transfer  ooeffiolents  for  an  aooeptable  flow- 
friction  power  expenditure)  (o)  reduced  thermal  stresses  within  the  fuel 
element)  euid  (d)  possibly  reduced  fabrication  costs  in  produotioa  quantities. 
Heat-transfer  and  flow-friction  oharacteristios  of  these  matrloes  are  of 
Interest  therefore,  for  future  reactor  developments  and,  moreover,  for  more 
oommonplaoe  applications  such  as  packed  absorption  and  extraction  columns) 
and  the  "peri odio- flow"  regenerative  heat  exchanger,  now  of  interest  in 
automotive  gas-turbine  application,  emd  (^. 

A  program  was  started  in  1949t  under  Office  of  Naval  Eesearoh.  sponsor¬ 
ship,  at  Stanford  University  with  the  general  objectives  of  investigating 
porous-body  heat  transfer  and  flow-friotion  behavior*  The  initial  studies 
were  conducted  with  woven-wlre  screen  matrices,  as  these  were  capable  of 
fairly  olear-cut  geometrical  definition  and  their  behavior  is  also  of  direct 
teohnioal  interest*  A  sphere-bed  matrix  also  was  included  because  sphere 
bods  have  served  as  a  "standard"  for  the  extensive  flow-friotion  inves¬ 
tigations  of  porous  media.  The  results  of  this  program,  throu£^  1932,  are 
available  in  references  (g),  (2),  and  (^a  A  more  recent  investigation  using 
packed  steel  wool  was  reported  by  Maroo  and  Han  (^» 

One  limitation  on  these  earlier  results,  from  the  point  of  view  of 
application  to  reactor  fuel-element  design,  was  that  they  were  limited  to 
a  Np  range  of  about  3  'to  1000.  At  the  suggestion  of  the  Eeactor  Engineering 
EivTslon  of  Argonne  National  Laborato:^,  cmd  with  flmanoial  support  from  the 
Atomic  Energy  Commission,  the  Office  of  Naval  Eesearoh  program  at  Stanford 
was  then  directed  to  increase  the  Eeynolds-number  range  up  to  100,000  for  the 
woven-soreen ^'geometry.  Consumatlon  of  this  objective  necessitated  a  complete¬ 
ly  new  experimental  teohnlque  and  test  setup*  The  fine-woven  screens  ranging 
from  3  to  60  mesh  with  wire  diameters  fxaa  7*6  to  41  mil  were  "modeled"  by 
orossed-rod  matrices  composed  of  0a373''iu~4iam  rods*  This  inorease  in 
geometrical  soa^e,  supplemented  by  increased  air-mass  veloolties,  provided 
the  desired  100-fold  inorease  in  N^* 

s 

Because  the  weaving  effect  of  the  screens  was  not  duplicated  with  the 
orossed-rod  matrices,  departures  from  oomplete  geometrioal  similarity  resulted* 
Either  porosi^,  p,  or  transverse  pitoh  x.,  oould  be  duplioated,  but  not  both* 
The  decision  was  to  duplioate  porosity)  am  as  a  consequenee,  small  oorreotions 
were  applied  to  the  crossed  rod-matrix  data  to  estimate  the  desired  woven- 
soreen  behavior  at  the  high  Eeqpnolds  ntanbers*  However,  heat- transfer  and  flow^ 
friotlon  behavior  of  the  orossed-rod  matrices  are  also  of  teohnioal  Interest* 
Consequently,  the  objectives  of  this  paper  are  to  report  the  behavior  of  both 
the  woven-soreen  and  the  orossed-rod  matrloes  in  the  %  ranges  from  3  'to  100,000, 
and  300  to  100,000,  respectively.  AdditlonsQ.ly,  it  is  proposed  to  describe  briefly 


1 

Underlined  numbers  in  parentheses  refer  to  the  Bibliography  at 
the  end  of  the  paper* 


the  different  test  techniques  used  to  achieve  the  low  and  the  high  test  ranges 
This  paper  constitutes  a  sxusBaiy  and  condensation  of  reference  (11) • 


RESULTS 

The  crossed<-rod  and  woven-screen  matrices  are  described  in  Table  1*  The 
dimensions  of  the  nominal  l6-mesh,  woven-screen  and  its  associated  crossed-rod- 
matrijc  model  are  shown  in  Fig.  1.  The  model  for  the  nominal  2U-^sh  woven^- 
screen  matrix^  described  in  Table  1,  is  shown  in  Fig.  ll;.  For  both  the  screen 
and  model  matrices  each  "screen  element"  was  orlsfnted  at  deg  with  its 
neighbor  so  as  to  simulate  a  "random  stacking''  in  contrast  to  a  "regular 
stacking,"  where  all  the  rod  elements,  for  instance,  would  be  either  parallel 
or  at  right  angles  with  each  other. 

Geometrical  Parameters 


The  geometrical  parameters  used  in  the  description  cf  either  a  wovexr* 
screen  or  a  crossed  rod  matrix  are  the  porosity  p,  the  flow  hydraulic  radius, 
rii,  the  rod  or  wire  diameter  d,  the  transverse  pitch  the  longitudinal  , 
pinch  ,  the  minimum  flow-area  ratio  cT,  and  the  heat-transfer  area  density,  a-*  « 
The  general  definition  of  the  hydraulic  radius  is 


rh  AfiojrA 

If  one  now  defines  an  average  flow  area, 

.  4 


‘flow 


P  A. 


frontal 


. (1) 

distinct  from  the  minimum  flow  area,  as 

. (2) 


since  by  definition 

. 

the  resulting  working  expression  for  evaluating  the  hydraulic  radius. is 

. 

It  is  to  be  noted  that  both  p  and  cc  can  be  measured  or  specified  without  ainblgulty* 


The  sketches  of  Fig.  1  describe  the  differences  between  the  crossed-rod 
and  wire-matrix  elements.  For  the  idealized  crossed-rod  geometry  shown,  the 
following  equations  may  be  derived  readily  from  Equations  (1)  to  (U)o 


p  ■  1  -  1P£/(U  x^) 
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(6) 
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(7b) 
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Note  that  porosity  is  a  fiuiction  of  only  for  the  crossed-^rod  matrix  from 
Equation  (5)  and  similarly  frmn  Equation  (7b) »  This  is  not  the  case 

for  the  wovKun^screen  matrix  because  of  wire  inclination  due  to  weavingo 
Equations  (7a)  and  (8)  apply  also  to  the  woven-soreen  matrixo 

Correlating  Parameters 

Two  Reynolds  numbers  will  be  used  for  the  purpose  of  correlating  the 
flow-friotion  and  heat“transfer  oharaoteriatioe 


®  G/vX,  and  Nj^^^  -  d 

It  follows  directly  from  the  definitions 


o  o  o  o  o  o 


®mai  “  ^^frontal^ 

together  with  Equations  ("ajTbjSj^)  that  generally  for  woven-wire  and  oroes- 
rod  matrices 


V%(d)  “ 

and  specifioaxly  for  the  crosaed-rod-matrix  geometry 
^E/*R(d)  *  ^  ^^t  “  »  • 


0  o  o  o  • 


.  .(10a) 


(10b) 


The  heat-transfer  parameter  which  will  be  used  is  the  oonventional 
Stanton  nvimber 

»3t  - 

Two  f^low^frlction  parameters  will  be  used,  the  Fanning  friction  factor 


o  o  o  o  o 


(11a) 


and  a  drag  coefficient 


1  2  2 
max/  P  i 

2g.  ^ 


(12a) 


For  essentially  constant  density  flow 


A 


frontal 


-6- 


Equations  (lla)  and  (l2a)  will  then  reduce  to  the  more  oonventloa8l*hydraullo 
form" 


In  Equation  (l2h)  Cj.  is  a  drag  ooeffioient^  or  pressure  ooeffiolenty  for  a 
single  soreen  of  thickness  g  «  2d:i^  . 

It  is  evident  that  oonversion  from  one  friction  parameter  to  the  other 
oan  be  made  from  the  general  relationship. 


(13a) 


and  for  the  orossed-rod  matrix  geometry 

%  ‘  2[j,r  (x,  . 

Equations  (lOb)  attd(l3b)  are  shown  graphioally  Im  Tig.  2. 

t 

As  mentioned  earlier,  the  small  geometrloal  dissimilarity  between  the 
woven-soreen  matrix  and  the  orossed-rod-matrlx  model  prompted  the  selection  of 
porosity  as  the  common  parameter.  The  observed  behavior  of  the  model  was  then 
"oorreoted**  to  estimate  the  behavior  of  the  woven-soreen  matrix.  The  oorreo- 
tlon  factors  are  described  in  the  following  text. 


From  the  flow  sketch  in  Fig.  1  it  is  evident  that  beoause  of  wire  inclina¬ 
tion  due  to  weaving,  the  boundary  layer  develops  over  a  longer  flow  length  than 
for  the  orossed-rod  element.  A  simple  boimdary-layer  analysis  indicates  that 
the  orossed-rod  element  would  have  hl^er  heat-transfer  ooeffioient  at  a  given 
mass  velocity,  beoause  of  the  shorter  perpheral  flow  length  for  boundary-layer 
development.  In  effect,  the  boimdary  layer  Reynolds  number  for  the  Inolined 
cylinder  is  greater  than  for  the  cylinder  normal  to  the  flow.  Thus,  for  the 
same  the  effective  Ep,  based  on  r.  ,  for  the  soreen  matrix  must  be  less  than 
for^^e  rod  matrlxi  and  analysis  indicated  that  cos  6  is  the  factor  to  apply. 


oos  B  E-  (orossed-rods)  •  E  (woven-soreen) 
—  R 


This  conclusion  was  oonflzmed  by  tests,  reported  in  (ll).  on  a  single 
cylinder  at  various  inclinations  up  to  34  deg.  As  Indioat^  in  Table  1,  the 
magnitudes  of  &  are  such  that  this  adjustment  on  Reynolds  nximber  is  less  than 
12  per  cent  for  five  of  the  six  matrines.  and  leas  than  2$  per  oent  for  the 
sixth.  Moreover,  slnoe  Eg.  oo  Ep  “0*37*  approximately,  even  the  25  per  oent 
adjustment  has  only  a  9  par  oent  influenoe  on  Eg^. 


-  7  - 


With  this  oorreotion  it  was  found  that  the  low  screen  heat-transfer 
behavior  correlated  very  well  with  the  hi^  orossed-rod  test  results* 

From  aerodynamic  oonsiderationf  and  also  from  a  oomparison  of  f/2  versus 
Wgt  Np  2/3^  the  oonolusion  was  reaohed  that  fora  drag  provides  the  major 
oontriDution  to  friction^  and  that  skin  friction  is  only  a  minor  factor*  For 
this  reason  a  oos  0  correction  is  not  warranted^  as  for  the  oase  of  heat  transfer. 
It  was  then  determined  hy  trial  that  the  direct  test  results  of  C-  versus 
with  porosity  as  a  third  parameter »  provided  a  very  good  correlation  of  the^l^ 
Hp/.\ crossed  rod-matrix  data  with  the  low  W./.x  woven-soreen-matrix  data. 

Bedause  the  model  x.  did  not  match  the  woven^aoreen  matrix  x.f  this  correlating 
prooedure«  in  effeoxy  resxilted  in  a  relative  shift  on  both  tne  friotion-faotor 
ordinate  and  the  Heynolds-number  absoissa  in  converting  the  orossed-rod  matrix 
Op  versus  H-/ to  f  versus  W-  for  the  woven-soreen-matrix  behavior*  The 
magnitudes  or  these  shifts  Is^lndloated  in  Fig*  2  by  the  ordinate  difference 
between  the  woven-soreen  matrix  of  specified  porosity  and  the  associated  orossed- 
rod  model  of  the  same  porosltyy  but  all^^tly  different  x^  It  may  be  noted  that 
the  friction  factor  and  Eeynolds-number  shifts  are  generally  small  except  for  the 
p  <■  0*602  matrix  where  they  are  quite  large* 

Heat-Transfer  and  Flow-Friotion  Results 


In  Figs.  3  to  8  are  presented  the  test  results  for  six  woven-soreen  matrioes 
ranging  in  porosity  from  83*2  to  60*2  per  oent*  The  correlating  factors  used  are, 
for  heat  transfer,  Wg.  H-  2/3  versus  Bp,  and  for  flmr-frlotion,  f  versus  Bp.  In 
each  oase  the  lower  B^  daxa  points  are  from  the  previously  reported  results^of 
Coppage,  (^*  The  hl^er  Bp  points  are  oorreoted  from  the  orossed-rod-matrix 
model  data,  as  described  previously*  For  heat  transfer,  the  Bp  was  multiplied 
by  008  0  to  acooimt  for  wire  inollnatlon  due  to  weaving*  The  xriotlon  factor 
f,  was  oaloulated  from  the  model  orossed-rod  matrix  0.,  and  Squatlon  (l3a),  using 
the  woven-wlre*4natrix  parameters  (r.t  0  >  P,  andff).  '‘'There  is  also  Involved  in 
this  -to-f  transformation  a  conversion  from  'bo  Bp,  using  Equation  (lOa), 

with  the  appropriate  woven-^re-matrlx  geometrloai  ^rameters* 

The  curves  representing  the  heat-transfer  behavior  (Figs*  3  to  8)  were 
oaloulated  from  the  emplrioal  equation* 

*Pr  »*E  . (W) 

where  Bp,  is  a  ’hnodifled"  Reynolds  number  as  follows 

f  -  0*96  for  B^  >1800 

f  -  1.155  -  0.0601  lotj^Q  BJ  for  B«  l800 
To  facilitate  the  use  of  Equation  (14)»  Big*  10  was  prepared,  giving 

K«/Bp  -(1-I^)/Fp  as  a  fiinotlon  of  Hp« 


s- 


The  six  heat-transfer  curves  from  3  ard  8  and  a  sever th  curve  also 

calculated  from  Squation  (14)  are  sximmarlzed  in  Fig«  9°  This  last  ourve^  for 
a  porosity  of  O039O9  provides  a  comparison  with  the  low  sphere  matrix  data 
reported  by  Coppage  and  London  (^o  The  agreement  is  surprisingly  good;  but 
somewhat  fortuitouso  because  of  the  great  difference  in  geometryo  Moreover; 
Denton* 8  results  for  fran  530  to  53;000;  which  can  bo  represented  by 


»3t  *Pr 


237  \ 


-0o31 


while  forming  quite  a  reasonable  extension  of  the  low  H-  sphere  data;  are  some 
56  per  cent  higher  than  the  curve  from  Equation  (I4)  "  50;000o  It  is 

recommended  that  Equation  (14)  be  restrloted  in  application  to  either  screen 
or  orossed-rod  matrices  In  the  range  of  porosities  of  55  *^0  ^5  por  oento  Mote 
that  the  application  to  randomly  stacked  crossed'^rod  matrices  will  be  moderately 
conservative  (by  less  than  10  per  cent  )  because  of  the  inoluslon  of  the  cos 
correction  factor  to  aoooxint  for  the  weaving  of  the  screen  matrloeso 


In  Flgo  11  are  presented  the  direct  test  restilts  of  C.  versus  for 
both  the  woven-soreen  and  crossed— rod  matrioeso  The  generally  good  agrenieat 
between  the  low  behavior  of  the  screens  and  the  hl^  ^B(d) 

determined  from  tnd  rod  matrices  supports  the  adequaoy  of  this^'type  of  correlation. 


An  empirical  equation  for  this  family  of  curves  was  determined  to  be 


Xo33.  w 

H(l) 


-0.33  .  SiSk 
P 


•  (15) 


The  adequaoy  of  Equation  (15)  in  representing  the  experimental  results  is 
demonstrated  in  Fig.  13 9  where  the  test  data  in  the  form  of 

0/  versus 

are  graphed.  The  agreement  with  the  line  from  Equation  (l?);  namely 


is  generally  within  I5  per  oent.  Also  on  this  graph  are  the  results  of  RomiS; 
et  al  (7) for  three  screen  matrices  and  these  data  also  support  the  oorrelatlon* 
Application  of  Equation  (13)9  as  for  the  heat-transfer  oorrslation  Equation  (14)9 
should  be  restricted  to  a  porosity  range  of  55  85  per  oento 

Because  of  possible  Interest  in  orossed-rod  geometries  for  reactor  fhel 
elements;  Fig.  12  is  yioluded.  Here  the  direct  test  results  are  presented  in 
the  form  of  Eg,  IL  2/3  f  versus  Ep«  It  is  recommended  that  this  representa¬ 
tion  bo  employod  xSr  porotitlos  la  the  ranM  60  to  83  per  cent  and  for  E-  from 
300  to  IOO9OOO  in  proforenoo  to  Equations  (I4)  and  (15)9  whloh  were  developed 
primarily  for  the  woven-soreen-type  matrix  and  therefore  inolude  the  small  cos  # 
oorreotlon  on  heat  transfer o 


J 
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EXPEBIUEKTAL  METHOD 


Th«  experimental  method  and  the  equipment  for  the  low  N-  oharaoteristios 
are  deeoribed  in  detail  by  Coppage  (^)o  These  tests  were  with  the  six  woven- 
soreen  matrices  described  in  Table  1  and  one  sphere  matrix,  porosity  >«  0o390, 
using  0a08l8--lno  lead  ballso  The  friotion  oharaoteristios  were  established 
by  oonventional  steady-flow^  isothermal  prooadures  and  the  heat  transfer  by  a 
transient  teohnlque  as  follower 

The  matrix,  initially  at  a  uniform  temperature,  is  heated  with  the  air 
flow  whioh,  at  "time  zero,’’  enters  at  a  constant,  hi^^er  temperatures  The 
temperature-time  history  of  the  fluid  leaving  the  matrix  is  recorded  during 
this  heating  process a  This  heating  ourve  is  oorpared  to  the  Anzelius  and 
Sohumtum  analysis  of  heat  transfer  to  an  idealized  perdue  body  and  the  convec¬ 
tive  heat-transfer  ooeffiolent  is  established  therebyo  The  oomparlson  with  the 
theoretical  solution  is  aooomplished  by  matohlng  the  maximiun  slope  of  the  heat¬ 
ing  ourve 0  The  reasons  for  seleotlng  this  procedure  in  lieu  of  other  methods 
of  oomparlson  are  discussed  in  reference  i^o 

The  high  Ip  oharaoteristios,  obtained  from  tests  on  orossed-rod  matrix 
models  of  the  woven-soreen  matrices,  were  determined  using  the  transient 
teohnlque  described  by  London,  Nottage  and  Boelter  (8) for  both  heat-amd  mass- 
transfer  investigations  and,  more  recently,  in  its  application  to  tube  banks 
(2)  and  (10) « 

Briefly,  the  teohnlque  consists  of  determining  the  time  rate  of  ooollng 
of  a  thermal  oapaoltor  of  the  appropriate  geometry  Immersed  in  the  air-flow 
streamo  In  Flgol4,  the  oapaoltor,  oonstruoted  of  copper  with  plastic  end  pieces, 
is  shown  as  one  rod  elmaent  of  the  matrixe  Allthe  other  rods  are  of  plastloe 
The  oapaoltor  oould  be  withdrawn  and  re-inserted  into  the  matrix  at  will,  with¬ 
out  interruption  of  the  airflow# 

Flgo  13  shows  the  section  of  duot  containing  the  test  matrix  along  with 
the  location  of  the  removable  thermal  oapaoltor#  With  the  airflow  established 
at  a  desired  rata,  the  oapaoltor  rod  was  withdrawn  from  the  matrix  and  heated 
in  a  small  eleotrio  oven  90  "feo  60  dog  P  above  the  duot  air-flow  temperature# 

The  thezmal-oapaoltor  rod  was  then  reinserted  into  the  test  matrix,  and  its 
tmnperature-time  history,  for  the  slower  oooling  rates,  was  determined  with  a 
Brown  indicator,  (25-ia  soalej  5^^1Iivolt  range?  l/2  seo  full-soale  response), 
noting  the  times  with  l/lO  seo  stop  watches  at  loO,  0e9t  O08,  0#?,  0e6,  and 
0a3  millivolts  above  the  temperature  of  the  reference  junction  in  the  air  stream# 
Thus,  the  test  ooollng  range  was  about  13  deg  F  starting  frexn  about  30  deg  F 
above  the  air- flow  temperature#  For  the  faster  ooollng  rates  a  Brown  Eleotronik 
recorder  was  used,  with  a  scale  setting  having  a  sensitivity  of  4  microvolts  and 
a  full-soale  response  of  better  than  1/2  seo# 


stream 


The  equation  for  the  temperature  of  the  thermal  oapaoltor  oooling  in  an  air 
,  with  a  pure  oonveotion  reslstanoe  at  the  surface  of  (l/hA)  hr  F/Btu,is 
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where  the  temperatures  are  referred  to  an  aiir— stream  datm  and  C,  the  themal 
capacity,  is  obtained  from  the  product  of  mass  and  specific  heat  of  the  copper® 
For  the  small  cooling  range  considered,  the  thermooouple  anf  response  is  linear 
with  temperature  so  that  the  minus  logarithm  of  the  emf  ratio  plotted  versus 
time  yields  a  straight  line  of  slope^/^^„  The  convective  heat-trsuiafer  coeffi¬ 
cient  is  calculated  from  this  slope® 

The  magnitude  of  C/A  for  the  thermal  capacitor  was  about  0o40  Btu/(F  sqft) 
(Table  l)o  As  the  range  of  measured  h  was  frcxn  3  to  100  Btu/(hr  sq  ft  deg  P)  the 
"time  constant"  for  the  thermal  resistanoe-capacitor  circuit  (c/hA),  ranged  from 
0®133  to  O0OO4  hr  (or  about  15  seo)®  From  this  last  magnitude  it  is  clear  why 
a  recording  potentiometer  was  required  for  the  faster  cooling  runs® 

For  the  crossed-rod  matrix  tests,  determinations  for  f  and  were  made 
from  Equations  (lib)  and  (I2b),  since  essentially  constant  density-flow  condi¬ 
tions  obtained®  Because  of  the  wide  flow  range,  it  was  necessary  to  measure 
the  pressure  drops  with  either  canpoxind  vertical-and-inolined  draft  gages,  or 
with  a  1-in-range  micromanometer  which  had  a  sensitivity  of  0o002  in®  of  water® 
The  measxired  pressure  drops  ranged  from  10  in®  down  to  OoOl  in.  of  water® 

As  a  preliminary  check  prior  to  all  these  tests,  pitot- tube  velocity 
traverses  were  made  upstream  from  the  test  section  in  order  to  demonstrate 
imiformlty  of  the  approach  flow,  and  also  provide  a  oheok  on  the  orifice 
flow  meters.  After  suitable  adjustment  of  the  trim  tabs  on  an  egg-crate  flow 
straightener,  located  upstream,  Z  1  per  cent  uniformity  over  96  per  oent  of 
the  flowsection  was  aohievedj  and  the  orifice  meter  also  was  checked  to  within 
1  per  cent®  This  flow  uniformity  oheok  was  made  at  only  one  air  rate,  but 
since  the  fxow  was  induced  from  the  laboratory  rown  throug^h.  a  smoothly  converg¬ 
ing  entrance  section,  and  the  duot  Heynolds  nvunber  was  always  in  excess  of 
30,000,  the  same  degree  of  xiniformity  can  be  reasonably  anticipated  for  the 
complete  test  range® 

The  estimated  experimental  uncertainty  intervals  associated  with  these 
two  transient  techniques  are  as  follows  ll) 


Uncertainty  Intervals 


Low  Np  tests,  a 

. .  ■  iL  per,  cgfll 

±15 

5 

2 

2 


Hig^  tests,  b 
_ a  per  cent 

i  3 

5 

2 

2 


a  Woven-screen  matrices 
b  Crossed-rod  matrices 
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From  the  tept  results  presented  graphloally  In  Figs.  3  to  8  and  Fig.  12p 
It  Is  olear  that  the  experimental  scatter  and  reproduoihllity  are  well  within 
the  foregoing  llmltso  Note,  however,  that  the  foregoing  estimates  of  exper¬ 
imental  uncertainty  do  not  Include  the  speculative  aspects  of  the  oorreotlons 
applied  to  the  oro8sed*°rod  matrix  behavior  to  predict  the  hi£^  I-  behavior  of 
the  woven-wlre  matrices.  These  Include  (a)  the  obliq,ue  flow,  oof  oorreotlon 
to  for  the  heat-transfer  data|  and  (b)  the  relative  shift  in  friction  factor 
f,  and  N„9  reported  In  Fig*  2,  as  was  suggested  by  the  versus  oorrela- 

tion  of  fig.  llo  ^ 

A  more  detailed  presentation  of  the  test  equipment,  the  teohnique,  and 
the  evaluation  of  experimental  uncertainties  is  oontalned  in  reference  (ll) . 


DISCUSSION 


In  all  the  tests  reported  in  this  paper  the  fluid  used  was  air  at 
approximately  room  temperature.  As  temperature  differences  for  heat  transfer 
were  always  small,  variations  of  properties  In  the  boundary  layer  were  neg¬ 
ligible*  In  application,  however,  some  consideration  must  usually  be  given 
to  the  influence  of  fluid  properties  variations* 

2/3 

A  Np  '  factor  Is  Included  for  the  Prandtl-nvuaber  effect,  even  though 
this  factor  was  essentially  constant  for  the  tests*  It  should  be  a  reasonably 
good  approximation  for  a  ^  range  of  0*6  to  3*  It  Is  not  applicable  for  the 
very  low  Frandtl  Number  or  liquid  metals.  Some  limited  experimental  results 
of  the  influence  of  Np^  in  crossflow  are  discussed  in  a  paper  by  Welsman  (12)* 


Two  Eeynolds  numbers  are  used  In  the  correlations,  Np  and  && 

defined  by  Equations  (9)«  The  first  of  these  mi^d^t  be  termed  an  equValent 
tube-flow  Eeynolds  number,  as  the  mass  velocity  Q  is  based  on  a  porosity 
average  flow  area,  and  r^^  oharaoterisos  a  flow-passage  dimension.  In  contrast, 
^B(d)  ^  parameter  oharaoteristlo  of  flow  over  a  body  of  dimension  d, 

expefienolng  an  oncoming  mass  velocity  0  *  From  the  point  of  view  of  the 

designer,  the  tube  flow,  or  Interior-flo^^oint  of  view  is  probably  the  most 
simple  to  apply  because  the  dimensions  r.  ,  p,  and  cb  are  more  closely  related 
to  the  design  problem*  In  contrast,  for^flow  friction  the  exterior  flow  or 
aerodynamio  point  of  view  is  most  useful  in  tying  together  the  orossed-rod 
and  woven-soreen  behavior.  For  design  oohvenlenoe,  after  the  matrix  geometry 
is  specified,  It  Is  reoammended  that  Op  versus  Sp/.\  be  converted  to  f  versus 
Np,  using  Equations  (lO)  and  (13) *  ^  ' 


It  Is  oomfortlng  that  the  f  oharaoteristlo  for  the  woven-soreen  geometxy 
at  the  high  Eeynolds  numbers,  as  predicted  from  the  model  tests,  line  up  so 
well  with  the  direct  tests  at  the  low  Np(Flg8*  3  to  8)*  This  agreement  is  not 
as  strong  a  verlfioation  of  the  C-  versus  oorrelation  as  desiredi  however, 

because  an  examination  of  Fig*  2  ^11  show  wai  only  the  p  ■  0*602  matrix 
Involves  a  strong  enough  relative  shift  in  transferring  to  the  f  versus  Np  basis. 
The  shifts  provided  for  the  p  ■  0*67?  and  0*725  matrloes,  while  noticeable,  are 
not  marked  enough  to  be  really  oonvinolng*  Moreover,  because  of  the  negative 
slope  of  the  f  versus  N-  oharaoteristlo,  these  Eeynolds  numbers  and  f  shifts  are 
partially  compensating*  In  effect,  while  the  0^  versus  oorrelation  of  both 


crossed— rod  and  woven-^wir©  maiiiicec  appear*  to  la  proa.i*i:'g9  furt^ier 

confirmation  for  porosities  in  the  ran,ge  of  5'  to  6?  per  cer*-  is  neededo 


I'he  correlating  (I4)  ar.d  (13}  fit  the  wove'-i-icreen-riatrix 

behavior  very  wellj  especially  when  the  ss-ide  range  of  fxoa:  5  "to  I'-OjODi 
is  taken  into  aocoanto  For  heat  transfer p  with  the  exception  cf  tne  range 


heloir  30  for  the  one  matrix  of  lo725  porcsitvj  agreevrent  is  wi 


“5  p0- 


oento  Neverthelessj  it  is  to  he  sti-ongly  empLar;i?ed  r-hat  while  these  equations 
are  adeq.uate  for  interpolating  withii';.  tt.e  test  poTO--’.il.y  rangep  they  carinct  ce 
safely  extrapolated  ontsida  of  the  range  'c55  P  <. -'£,85.  In  spite  of  thii-: 
cautionp  heoaase  of  the  lack  of  basic  design  dstap  tee  engineer  aay  be  forced 
to  make  extrapolations  not  only  to  porosities  ontside  this  rangej  but  also  to 
other  geometries*  Some  support  for  this  speculative  prooednrep  for  the  heat- 
transfer  behavior p  is  provided  by  the  fact  that  the  Kpher*  behavior  on  Pig.  9 
is  predicted  q_ulte  well  by  Eq-uation  (14)  for  h?.000  and  even  at  a  - 
50,000  the  extrapolation  provided  by  Equation  fid)  is  only  5^  per  cent'  below 


the  resvilts  of  Denton  (_^( 
trai 

14). 


_  Comparisons  of  the  low  woven-acreen  heat 
transfer  are  in  approximate  agreement  with  that  of  bed?,  of  granular  materials 


The  heat-transfer  behavior  for  flow  no.mal  to  tube  banks  i^s  for  the 
two  staggered  arrangements  with  a  longitudinal  pltoh  of  unityp  as  for  the 
orossed-rod  matrices p  are  also  in  fair  agr'sement  (21  to  30  per  cent  lower) 
with  Equation  (14)«  The  transverse  pitch  for  these  geometries  are  x^,  -<  2.00 
and  1.50  corresponding  to  porositiesp  fran  Equation  115} j  of  vo607  and  0»477> 
respeotively.  Further  oompariaons  of  this  nairare  are  reported  by  Weiaman  (12} » 

The  behavior  reported  here  for  oroseed-xod  and  woven-soreen  matrices  are 
for  ’'random  stacking"  in  contrast  to  "reg'ular  stackingp"  where  the  rods?  ox- 
wires  are  either  parallel  or  at  right  angles  with  each  other.  For  fuelelements, 
the  "regular  stacking"  may  prove  to  be  more  attractive  from  the  viewpoint  of 
fabrloation.  Moreovexp  there  is  some  evidence  to  believe  that  the  heat- 
transfer  behavior  may  be  improved  by  about  I'J  per  cent  with  the  sam*  or  lower 
frloti(r::ip  for  the  "staggered  regular  slacking"  where  flow  self-ahielding  is 
not  as  pronounced  as  in  either  the  "in-lfnce  regular  stacking"  or  "rendosE  stack¬ 
ing"  geometries.  These  questions  are  currently  nr.ii;:-  investigation. 


COECLDSIOirS 

The  following  oonolusions  result  from  this  investi^tioas 

1  Adequate  basic  heat- transfer  and  flow-friotion  data  are  now 
available  for  the"rardom  stacked"  woven-soreen  and  orossed-rod  types  of  matrices 
over  a  wide  range  of  Eeynolds  numbers  and  a  porosity  range  from  55  fo  85  per 
cent. 


2  The  fact  that  two  markedly  different  transient  test  techniques 
were  employed  to  cover  the  Reynolds  number  rangep  and  that  substantial  agreement 
results  in  the  overlapping  regionp  provides  confidence  in  the  validity  of  the 
results. 
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3  Extrapolation  of  these  results  to  the  prediction  of  the  heat- 
transfer  and  friotlon  hehavlor  of  other  porovis  body  geometries  la  ri8lgr«  How- 
evBTf  for  the  heat- transfer  behavior  at  lease,  some  limited  support  for  this 
procedure  is  in  evidence o 

4  The  transient  technique  used  for  the  orossed-rod  matrices  has 
proved  to  be  both  simple  and  accurate <>  Wider  applicatlcn  to  other  convective 
heat-transfer  Investigation  is  suggestedo 

5  For  the  orossed-rod  matrix  geometry,  remaining  investigations 
include  the  '"regular  stacked"  arrays  of  both  the  "in-line"  and "staggered" 
variety  and  a  wider  range  of  porosities  for  both  the  "random"  and  "regular 
stacked"  gemnetrleso 


ACOOlLElXaSIEENTS 

The  work  reported  here  was  accomplished  as  part  of  a  program  at  Stanford 
University  sponsored  by  the  USE  Office  of  Naval  Research,  the  Bureau  of  Ships, 
the  Bureau  of  Aeronautics 9  and  the  ASOo 

The  authors  are  particularly  appreciative  for  the  teohnloal  assistance 
provided  by  Profo  Wo  M«  Kays  and  to  Kto  Ao  imorosl,  formerly  at  Argonne 
National  Laooratory;  for  suggesting  the  prograao 

Stanford  engineering  students  Mto  Ao  Ro  Ray  and  Ifr.  M,  So  Shihabl 
assisted  in  the  test  work,  oaloulatious,  and  draftingo 

The  Harrison  Radiator  Division  of  General  Motors,  throu^  Mr.  John  Godfrey, 
provided  the  very  excellent  mioromanometer  used  for  the  friction  factor  determina¬ 
tions  at  very  low  pressure  drops o 


BIBLIOORAPHT 


1  "The  Periodic-Flow  Regenerator  -  A  Summary  of  Design  Theory"  by 
Jo  B.  Coppage  and  A.  Lo  London,  Trans.  ASMS,  vol  75,  1953,  PP«  779*'787« 

2  "Compact  Heat  Exchangers,"  by  W.  M.  Kays  and  A.  L.  London, 

National  Press,  Palo  Alto,  Calif o,  1955;  PPo  19~20,  144~148* 

3  "Heat  Transfer  and  Flow  Friction  Charaoteristios  of  Porous  Media," 
by  Jo  Ee  Coppage,  PhD  Thesis,  Stemford  University  (1952).  Also  issued  as 
Stanford  University  Teohnloal  Report  No«  16,  Office  of  Naval  Research  (NR-O9O- 
104),  December,  1952o 

4  "Heat  Transfer  amd  Flow  Friction  Charaoteristios  of  Porous  Media," 
by  Jo  Bo  Coppage  and  Ao  Lo  London,  Chemical  Engineering  Progress,  Feb.,  1956, 

PPo  57°*F  to  63-Fo  This  is  a  summary  of  reference  (3). 


-14- 


3  "An  Investigation  of  Conveotlon  Heat  Transfer  In  a  Porous 
Medliasy"  by  8.  M.  Maroo  and  Lo  3.  Han,  ASME  Paper  Ho.  35~-d-104* 

6  "The  Heat  Transfer  and  Pressure  Loss  in  Fluid  Flow  throxi^ 
fiandottly  Paokod  Spheres |"  hy  W.  H.  Denton,  C.  H.  Robinson,  and  R.  S.  Tibbs, 
HPO-35  June  28,  1949*  Available  throu^  ABC  Teoh.  Info.  Div.  ORHL. 

7  "Heat  Transfer  and  Pressure  Drop  Charaoterlstlos  of  Four 
Regenerative  Heat  Exohanger  Matrioos,"  by  F.  B.  Romle,  et  al,  ASMB  Paper 
Ho.  5I-SA-34. 

8  "Determination  of  Unit  Conduotanoes  for  Heat  and  Mass  Transfer 
by  the  Transient  Method,"  by  A.  L.  London,  H.  B.  Hottage,  and  L.  M.  E.  Boelter, 
Industrial  and  Engineering  Chemistry,  vol.  33f  April,  1941»  PP*  467-473* 


9  "Heat  Transfer  and  Frlotlon  Charaoterlstlos  for  Qas  Flow  Hoxmal 
to  Tube  Banks  -  Use  of  a  Transient  Test  Teohnlquo,"  by  W.  M.  Hays,  A.  L.  London, 
and  R.  K.  Lo,  Trans.  ASMB,  vol.  76,  1934»  PP*  387-396. 

10  "Average  and  Local  Heat  Transfer  for  Crossflow  through  a  Tube 
Bank,"  by  R.  A.  DeBortoll,  R.  E.  Orimble,  and  J.  B.  Zerbe,  ASME  Paper  No* 

55-SA-5I. 

11  "Heat-Transfer  and  Frlotlon  Charaoterlstlos  of  Soresn  Matrloes 
at  Hi^  Reynolds  Numbers,"  by  Long  Sun  Tong,  PhD  Thesis,  Stanford  University 
(1936;.  Also  issued  as  Stanford  University  Teohnloal  Report  No.  28,  Offloe  of 
Naval  Research  (IR-O63-IO4) 9  April  1936* 

12  "Effect  of  Told  Volume  and  Prandtl  Modxdus  on  Heat  Transfer  in 
Tube  Beuiks  and  Packed  Beds,"  by  J.  Welsman,  AIChB  Journal,  vol  1,  September, 
1935,  P*  342. 


•15- 


Flgo  1 

Flgo  2 


Pig  3 

T±go  4 
Pigo  5 
Pig«  € 

Fig.  7 
Pig*  8 
Pig.  9 


EEAT-TRAirSFER  AND  PLOW~FBICTION  CSABACTERISTICS 
OF  WOVSH-SCRESN  AND  CBOSSBD-ROD  MATRICES 


CAPTIONS  FOR  ILLUSTRATIONS 


Geranetry  of  woven'-soreea  and  orosaed-rod  matrioas  (Dl  •mansions 
are  for  the  nominal  l6<HBesh  soreen  matrix  and  oorre^  ing 
orossed-rod  model*  See  Table  1*  Note  that  worea-sc  matrices 
may  have  either  greater  or  leas  than  imity  as  a  result  of 
orimping  or  flattening  respeotively.) 

Ploe-friotion  factor  and  Reynolds  nximber  ratios  for 
orosaed-rod  matrioas 
Prom  Bquatsioas  (l3b)  and  (lOb). 

+  denotes  orosaed~>rod  model  conditions 

•s,  0  denotes  woven-wire  matrix  oonditlonB  oaloulated  from  Equations 
(l3a)  and  (lOa)* 

(Ordinate  differences ,  for  a  woven-^re  matrix  and  its  oorresponding 
modelp  indicate  the  relative  shift  Involved  la  converting  from  the 
0.  versus  N-/^\ correlation  to  predict  the  f  versus  N^  behavior  of 
screens  at  ni^  N^  shcrwn  in  Pigs»  3  to  8*) 

Woven-soreen  matrix  oharaoterlstlos  N-.  N-  and  f  versus  N 

(See  Table  1  for  details  on  the  geometry*^  roruslty  p  •  0*832. ) 

Woven-soreen  matrix  oharaoterlstlos  Eg.  N_  and  f  versus  Np 

(See  Table  1  for  details  on  the  geomet^«^Foroiiilty  p  0*817.)^ 

Woven-soreen  matrix  oharaoterlstlos  Ng.  Np  and  f  versus  N- 

(See  Table  1  for  details  on  the  geometry*  ^Porosity  p  ■  0*766*7 

Woven-soreen  matrix  oharaoteristlos  Ng^  Np^  and  f  versus  N|j 
(See  Table  1  for  details  on  the  geometry*  Porosity  p  ■  0*72$*) 

Woven-soreen  matrix  oharaoterlstlos  Eg.  Np  and  f  versus  Np 
(See  Table  1  for  details  on  the  geomet^*  ^Poruslty  p  ■  0*675*7 

Woven-soreen  matrix  oharaoterletlos  Eg.  ^  versus  Np 

(See  Table  1  for  details  on  the  geomet^*  ^Porosity  p  ■  0*602*)" 

Summary  of  heat-transfer  oharaoteristlos  for  woven-soreen  natrioes 
*St  *Pr  h 

(Lines  are  oaloulated  from  correlating  Equation  (14) «  These 
are  the  same  as  the  lines  shown  in  Figs*  3  to  8  for  the 
soreen  matrioes*  Alee  shown  is  a  line  oaloulated  for  a 
porosity  of  0*390  (oorresponds  to  the  porosity  of  a  random 
packed  sphere  matrix  tested  by  Coppage  *) 
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*E*/*E.  ■  “  Fp)/Pp  vePTOB  I« 

Supplaaentary  graph  to  uaa  with  Bciuation  (I4) 


Priotion  oorralation  for  hoth  orossad-rod  and  wovan-aoraan 
■atrioaa  0^^  veraua  *5(4) 

and  f  Taraua 


Orosaad-rod  matrix  oharaotariatloa 


2/3 


Taat  of  frlotlon  oorralatlng  Bqnatlon  (13) 

Oroaaad-rad  matrix 

(Modal  for  p  ■  Oo725f  24  maah  wowan-^ra  matrix  daaorlbad  in 
Ti^la  la  Modal  matrix  oonaiata  of  5  "aoraans"  or  10  rowa  of 
roda.  The  thazmal  oapaoitor  alaaiant  ia  looatad  in  tha  front 
raw  of  tha  fourth  "aoraaa”  in  tho  flow  dlrootlona) 

Oroa»-rod  matrix  la  taat  duot 

Spaolfioatlona  for  Larga  Soala  Modal  and  Soraan  Matrloaa 


